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We present results from the first fully nonlinear numerical calculations of the head-on collision of 
two unequal mass black holes. Selected waveforms of the most dominant ^ = 2, 3 and 4 quasinormal 
modes are shown, as are the total radiated energies and recoil velocities for a range of mass ratios 
and initial separations. Our results validate the close and distant separation limit perturbation 
studies, and suggest that the head-on collision scenario is not likely to produce an astrophysically 
significant recoil effect. 
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The collision of two black holes is expected to be one of the most interesting and important astrophysical sources 
of gravitational radiation. The strong bursts of radiation emitted from the interaction should be detectable by the 
next generation gravitational wave observatories, and will provide explicit evidence for the existence of black holes, 
assuming the characteristic signatures of the holes are observed. To compliment the anticipated observations, several 
calculations of binary black hole interactions have been performed in recent years, inspired in part by a broader 
Grand Challenge effort to supply waveform templates that are essential for the analysis of observational data. These 
calculations have been based on axisymmetric numerical relativity computations of the head-on collision of two initially 
time-symmetric and boosted black holes as well as perturbation studies which confirm and elucidate the 

numerical results in the close limit approximation (CLAP). However, the fully nonlinear numerical studies performed 
until now have been restricted to the simplified case of equal mass black hole systems. To extend these calculations, 
we have recently introduced a new class of body-fitting coordinate systems (class I in the notation of Ref. ) that 
overcomes the coordinate singularity problems in previous work, fn this paper we utilize these new coordinates to 
carry out the first numerical studies of axisymmetric head-on collisions of unequal mass black holes. In addition to 
expanding the catalog of computed energies and waveforms, numerical investigations of these systems can address 
whether gravitational radiation from the collisions is likely to carry an astrophysically significant component of linear 
momentum, and whether substantial recoil velocities can be imparted to the final black hole. 

The initial data used in our simulations are time-symmetric analytic solutions that generalize the equal mass Misner 
data adopted in previous numerical calculations Applying the isometry operation systematically across each 
of the throats for the infinite image poles, the conformal factor can be written in cylindrical coordinates as 
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where we have defined /i„,i = /in, 2 — ^i^o for even n, ^n.i = fJ-i + (n — l)/io and fin.2 — fJ'2 + [n — l)/^o for odd n, and 
fiQ = {ill +/i2)/2. The initial data is thus uniquely specified by two parameters, fii and fi2, which determine the black 
hole masses and separations, analogous to the single Misner parameter in the equal mass case. The physical content 
of the initial data is characterized by the total (ADM) mass of the spacetime derived as the asymptotic leading order 
term in Madm = — 2r^9.r* 
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the proper distance separating the black hole throats along the [p = 0) axis 
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and the bare masses of the two holes as suggested by Lindquist ||J^ 
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The numerical solutions are parameterized by the ratio of bare masses Mr = mi/m2 with the convention mi < 7712, 
and the proper separation between the throats L/M normalized to the maximum of the two masses. Here we simply 
define M = Madm/{^ + Mr) for consistency with previous work. The code is applied to collide black holes with 
initial separations up to ~ lOM with mass ratios Mr = 1, 0.75, 0.5 and 0.25. These choices are determined mainly by 
stability issues, which require (/ii, /i2) ^ 3.0 as the upper bound 0|, and accuracy issues which require (/ii, ^12) ^ 0.5 
to resolve the very weak higher order {£ > 3) multipole waveforms, although an even greater range of parameters 
can be evolved accurately in the stronger £ = 2 signals. For each of the cases presented, the code is stable enough 
to maintain accurate solutions to well past IOOM^^jm, which is more than enough time to monitor the complete 
wave signals with current parameters. Although results are presented for Mji > 0.25, we expect that this range of 
parameters will produce maximal or near maximal recoil velocities, a conclusion supported by the CLAP |0] and 
Fitchctt 1^ calculations which predict maximum recoil at 0.25 and 0.38 respectively. 

The method we use to compute waveforms is based on the technique developed by Abrahams and Evans |^ (and 
applied in Rcf. to extract the Regge- Wheeler perturbation functions and to construct the gauge invariant Zerilli 
function ip. Fig. ^ shows a characteristic sample of generated waveforms for the particular case of /ii — 2.5 and 
/i2 = 2 corresponding to Mn — 0.66 and L/M — 7.45. Also shown are the quasinormal mode fits from perturbation 
theory using a linear combination of the fundamental and first overtone for the final black hole of mass Madm- The 
fits nicely match both the frequency and damping rate in the late-stage ringing phase. In Fig. ^ we show the total 
energy radiated in the dominant £ = 2 mode as a function of initial separation distance for different mass ratios. 
For comparison, we note that the perturbation results of Davis et al. (DRPP) when rescaled by the reduced 
mass, yields E/Madm ~ 0-0104^^/Af^^j^^ « 0.0104Af|,/(l + Mr)^ for the radiated energy in the large separation 
limit, where is the reduced mass. Our numerical results are in remarkable agreement with this prediction, which 
suggests that the radiated energy (when normalized to the ADM mass) decreases as and yields maximum energy 
emissions of (3, 7, 10 and 13)xl0~'' for the Mr = (0.25, 0.5, 0.75 and 1.0) cases respectively. In the opposite close 
separation limit, our results are also in excellent agreement with the CLAP calculations of Andrade and Price [Q, 
and the particle limit {Mr — > 0) collisions from finite initial separations based on the Green's function integration by 
Lousto and Price |l2| . This last curve appropriately bounds the sequence of data points from our numerical relativity 
calculations. 

As indicated by Fig. |^, the quasinormal modes are strongly excited, and the black hole oscillations are well described 
by essentially the quasinormal modes (and particularly the £ = 2 mode) of the final black hole despite the fact that 
the systems considered are not small perturbations. In this regard, the collision of unequal mass black holes is 
similar to that of equal mass systems, except now the breaking of equatorial symmetry allows for odd mode multipole 
components. The mixing of adjacent multipole modes gives rise to a nonvanishing flux of linear momentum along the 
z-axis which can be evaluated using the Landau-Lifshitz pseudotensor p3| , p^ 
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As an additional check we have also computed the momentum flux from consecutive Zerilli components 



dt IGiT f^\j {2£+l){2£ + 3) dt dt ' ^' 

and we find generally good agreement (typically 10 - 30% relative differences) between the two calculations. As 
a result of the momentum emission from odd/even mode mixing, the final black hole will acquire a recoil velocity 
= —M^^j^j J {dP^ /dt)dt, opposite in direction to the momentum fiux of the waves. 

The radiation of momentum by gravitational waves is of interest in work relating to active galactic nuclei, quasars 
and even archetypical galaxies since highly dynamic black holes and coalescing binary black hole systems may abound 
in galactic discs as well as in the centers of galactic nuclei. Because the efficiency of momentum radiation emission 
is not known precisely, the stability of these systems remains in question due to the generation of a recoil velocity in 
the affected black holes, which, for sufficiently asymmetric configurations, may be large enough to break black holes 
free from the host galaxy. If so, gravitational radiation effects will have considerable observable consequences for 
astrophysics and cosmology (such as the depletion of black holes from host galaxies, the disruption of active galactic 
core energetics, and the introduction of black holes and stellar material into the IGM) since the processes generating 
the momentum emission, i.e., stellar core collapse and binary mergers, are likely to be fairly common. 

Although fully nonlinear numerical simulations of asymmetric black hole systems have not been performed to 
now, several quasi-Newtonian and perturbation calculations have been carried out over the past two decades to 
estimate the magnitude of this recoil effect. For example, Bekenstein derived an upper bound of Vr < 300 
km/sec for the recoil velocity from nonspherical stellar core collapse; Moncrief found Vr = 25 km/sec for small 
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non-spherical distortions in Oppenheimer-Snyder collapse models; and Nakamura and Haugan (NH) computed 
Vr = \Ap\/M ~ 262(TOi/Af)^ km/sec for a test particle of mass mi plunging from rest at infinity into a Schwarzschild 
black hole of mass M, with mi -C M. Extrapolating the NH result to the comparable mass limit (mi ^ M), we 
replace the mass dependency with Fitchett's 1^ scaling formula, (mi/M)^ /(Afi?) — M]f {l — Mji){l + Mji)~^, which 
has the desired perturbation and equal mass limits. Considering the maximal value oi f — 0.01789 at = 0.38, the 
perturbation result predicts a maximum recoil velocity of about 5 km/sec for the head-on collision scenario. 

However, numerical calculations are needed to adequately resolve the equal or near-equal mass cases, and to 
properly model nonlinear effects from internal dynamics such as tidal heating, and the re-absorption and beaming 
of gravitational waves. We plot in Fig. || the absolute value of the recoil velocities found in our simulations by 
evaluating the integrated momentum flux from the Landau-Lifshitz pseudotensor across a 2-surface of radius 
15Madm- Despite the fact that we are unable to accurately evolve black hole collisions with initial separations 
greater than about lOM, Fig. |3| does suggest a maximum recoil on the order of ~ 10 — 20 km/sec, in rough 
agreement with the various perturbation and quasi-Newtonian estimates. At sufficiently small initial separations, 
our results also agree nicely with CLAP, especially when considering the smallness of this effect and the estimates of 
numerical errors based on resolution studies (less than 10% variation between fixed detectors in the 200 and 300 radial 
zone evolutions) and extractions of the momentum across differently positioned detectors (roughly 10-30% variation 
between detectors at radii 15 and 25 Madm)- In comparison, the escape velocity from galactic structures with radius 
R and total mass M is of order yjGM/R, which can vary from anywhere between about 200 km/sec for the less 
massive compact dwarf galaxies to about 1000 km/sec for giant ellipticals (though black holes in the galactic disk and 
far from the core would require substantially less recoil in the direction of galactic rotation to achieve escape velocity) . 
We thus conclude that the head-on collision of black holes is not likely to produce an astrophysically significant recoil 
effect. 

We wish to thank Z. Andrade, R. Matzner, R. Price and L. Smarr for many helpful discussions, and especially ZA 
and RP for also making available their CLAP data for comparison. PA is grateful to the Albert Einstein Institute for 
their hospitality during which part of this work was conducted. The calculations were performed on the C90 at PSC 
and the Origin2000 at AEL 
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FIG. 1. 

Zerilli waveforms for the /ii = 2.5 and 112 — 2 case, corresponding to a mass ratio = mi/m2 = 0.66 and proper 
initial separation L/M — 7.45. The ^ = 2, 3, and 4 modes are displayed normalized to the ADM mass, and extracted 
at a radius of 2QMa£)m- Also shown with dotted lines are fits to the fundamental and first overtone of the quasinormal 
modes from perturbation theory. The simulations are run at a grid resolution of 200x70 (radial x angular) zones. 
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FIG. 2. 

The total energy radiated in the dominant (. = 2 mode from the head-on collision of two black holes as a function 
of separation distance L/M between the holes, where M = MadmI{^ + Mb) is the approximate mass of the larger 
black hole, and — mi/m2 is the ratio of bare masses. The energies are computed as the integrated flux across a 
2-surface of radius IhMADM, and rescaled by the ADM mass and the perturbative (toi/to2 ^ 1) behavior. Four level 
curves are shown for different ratios of bare masses Mn = 1, 0.75, 0.5 and 0.25, together with the corresponding CLAP 
results from Andrade and Price (AP) the reduced mass extension of the Davis et al. (DRPP) perturbation 
calculation, and the Green's function integrations of Lousto and Price (LP) [T^ . 
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FIG. 3. 

The absolute value of the recoil velocity is plotted together with the CLAP results from Andrade and Price (AP) 
in units of km/sec for the three unequal mass level curves (Af_R — 0.75, 0.5 and 0.25). Due to the difficulty in 
resolving the weaker high-order multipole components (i > 3) in the close separation limit, the numerical data is 
shown only for (/ii,/i2) ^ 0.5, or roughly L/M ^ 2. The results here (and in Fig. ||) are from simulations run at a 
grid resolution of 300x70 (radial x angular) zones. 
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